Two-dimensional ion trap lattice on a microchip for quantum simulation 
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Using a controllable quantum system, such as trapped ions, it is possible to simulate highly com- 
plex quantum systems efficiently, overcoming an in-principle limitation of classical computing. So 
far, no dedicated architectures for the simulation of two-dimensional spin lattices using ions in radio- 
frequency ion traps have been produced, limiting the possibility of carrying out quantum simulations 
on a large scale. We report the operation of such an architecture integrated in a microchip. Our 
scalable microfabricated architecture allows for trapping of ion lattices with coupling strengths be- 
tween neighbouring ions sufficient to provide a powerful platform for the implementation of quantum 
simulations. To realise this device we developed a specialist fabrication process that allows for the 
application of very large voltages. We fabricate a chip containing an array of closely spaced micro- 
traps, demonstrate reliable trapping of ytterbium ion lattices, measure long ion lifetimes, perform 
rudimentary shuttling between lattice sites and show the suitability of the architecture for quantum 
simulations. We are also able to deterministically introduce defects into the lattice, adding further 
classes of quantum simulation that can be performed and allow for simulations with spins greater 
than 1/2. Another step-changing application consists of the generation of two-dimensional cluster 
states. 



An analogue quantum simulator is a device in 
which the Hamiltonian of a complicated many-body 
system can be realised and its properties measured. 
Such a device would enable investigation of quantum 
phenomena currently beyond that of classical com- 
putation or direct measurement [1, 2]. The appli- 
cations for a quantum simulator are numerous and 
far reaching, including, among others, the simula- 
tion of quantum magnetism, high-temperature su- 
perconductivity, the fractional quantum Hall effect 
and synthetic gauge fields [4-7]. Quantum simula- 
tions may revolutionise our understanding of many- 
body systems. 

Proof of principle analogue quantum simulators 
have already been demonstrated using ultra-cold 
quantum gases [8], photons [9], trapped ions [1, 2] 
and NMR [10]. Of the various systems trapped ions 
are one of the most promising due to the demonstra- 
tion of unparalleled quantum control and isolation 
from the environment [11]. 

While some seminal small scale quantum simu- 
lations have been performed within radio-frequency 
(rf) ion traps [12-18], they suffer limitations when 
trying to scale to a larger number of particles. The 
harmonic potential along the weak axis can result 
in varying ion-ion separations and coupling rates. 
Furthermore, ions at the rf nil naturally form a 1- 
dimensional string making the simulation of arbi- 
trary lattice geometries experimentally demanding 
[19]. A 2D Ising-type interaction has been observed 
with Penning traps [20], but the rotating crystal 
makes individual ion addressing and readout exper- 
imentally challenging, and the lattice geometry is 



limited to the naturally forming Wigner crystal. 

A different approach constitutes the proposal to 
use a two-dimensional lattice of rf ion traps [21-23]. 
Each lattice site contains a single ion which inter- 
acts with neighbouring ions via the Coulomb force. 
To date, experimental progress towards such a lat- 
tice has been limited to trapping dust particles and 
ion clouds above PCB boards [24] or wire meshes 
[25]. Optimised two-dimensional lattice geometries 
for the implementation of quantum simulations have 
been investigated focusing either on the topology of 
a single rf electrode to maximise the ratio of the 
ion-ion coupling rate to the decoherence rate [26] or 
focusing on maximising the electric field curvatures 
of individual trapping sites which can contain non- 
intuitive isolated rf and static voltage electrodes [27]. 
Here, we report on the microfabrication of a lattice 
of ion traps on a microchip and demonstrate trap- 
ping of a lattice of single Yb+ ions, the deterministic 
trapping of multiple ions at lattice sites and rudi- 
mentary ion shuttling between lattice sites. Addi- 
tionally this microchip offers a suitable architecture 
for 2D cluster state generation [3] and a possible av- 
enue for the development of quantum information 
processing [21]. As our device is an ion trap chip 
[28] it is scalable to much larger ion lattices. 

The layered structure of a commercial silicon-on- 
insulator (SOI) wafer was used by Britton et al. 
[29] to fabricate a linear ion trap permitting the ap- 
plication of voltages of up to 150 V. We have de- 
veloped a different fabrication process using a SOI 
wafer to fabricate a planar ion trap with a signifi- 
cantly increased trap depth whilst maintaining large 



ion-electrode separations using a combination of re- 
cessed electrodes and the application of voltages in 
excess of 1 kV. We describe both of these fabrication 
advances in more detail below. 

We utilise the layered SOI wafer to form a planar 
trap electrode structure, with a 2D electrode pat- 
tern but instead of the ground electrodes lying on 
the same plane as the rf electrodes they are recessed 
from the chip surface. This recess was achieved by 
using the SOI handle layer as rf ground, shown in 
Fig. lA and B. By using recessed electrodes, as 
shown in Fig. IB, the trap depth for a given ion- 
electrode separation and voltage can be substantially 
increased. 

In operating a microfabricated ion trap, restric- 
tions on the rf voltage that can be applied limit trap 
depth and require small ion-electrode distances. Low 
fiashover voltages have been a problem in previously 
reported microfabricated ion traps [28, 30], limiting 
obtainable ion life times and secular frequencies. To 
ensure efficient operation of our lattice ion trap we 
have developed a fabrication process that allows for 
high voltages, greatly reducing constraints on the 
trap design. We chose a thick oxide layer thus in- 
creasing the path length between the electrodes and 
the handle layer as well as reducing capacitance be- 
tween the rf electrode and ground. However this 
would only result in a modest increase in breakdown 
voltage. We were able to obtain a much larger in- 
crease in breakdown voltage by using a specialized 
fabrication process. In order to explain this pro- 
cess it is important to note that voltage breakdown 
will usually occur via insulator surfaces connecting 
two conductors rather than through the insulator 
bulk. We selected a particular SOI structure that 
was fabricated by wafer bonding two wafers with 
5 /im thick oxide surfaces to form a 10 /im thick 
buried oxide layer. Using a buffered HF etch, last- 
ing 130 — 140 minutes, the oxide layer was etched 
to expose the handle layer. As there is an increased 
etch rate on the interface where the two oxide lay- 
ers were bonded together a highly anisotropic etch 
results laterally under the electrodes. Rather than 
obtaining the usual approximately straight etch pro- 
file, we were able to obtain the V-shaped under- 
cut as shown in Figs. IC and 2B. This V-shaped 
etch profile substantially increases the path length 
between the electrodes and ground resulting in ex- 
tremely high breakdown voltages as it increases the 
effective distance for surface fiashover from 10 /im 
to up to 120 /im. Electrical fiashover measurements 
were performed on SOI test samples to determine 
the voltage that can be applied (see Methods). The 
mean fiashover voltages were measured to be \dc = 
1298(5) V and V^/ = 1061(32) V. These values are 
one or two orders of magnitude higher than in pre- 
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Figure 1: (A) The electrode design for the 
microchip. The rf electrode is shown in red, 
compensation electrodes are shown in gold and the 
ground electrode is shown in green. The hexagon 
geometry is shown in the insets, with a polygon 
radius r = 125 /im and the distance from the outer 
polygon centre to the rf edge, labeled w, is 410 /im. 
(B) Cross-section drawing of the trap geometry, 
showing the SOI layered geometry. The close-up 
inset shows the deep V undercut achievable with 
our fabrication process. (C) Two schematics of the 
etch profile, a typical isotropic etch is shown on the 
left, showing a slight undercut into the buried 
oxide. The right picture shows the profile of our 
optimised etch. The preferential etch rate along the 

bond face results in a deep V-shaped undercut. 
The ion height relative to the trap electrodes is also 
shown. 



viously fabricated microfabricated ion traps [28, 29]. 
In fact these voltages correspond to typical values 
applied in traditional ion traps made using large 
metallic rods illustrating the impressive performance 
of our microchip. 

The microchip consists of 29 traps arranged in a 
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triangular lattice, with each ion having up to six 
nearest neighbours with an ion-ion separation of 
270.5 jam. The ion-electrode separation is 156 jam. 
Due to the recessed ground electrodes the ion height 
from the top of the trap surface is 116 /im, as seen 
in Fig. IC. 




Figure 2: (A) An SEM of a finished microchip. 
The inset shows a close up of two of the hexagonal 
traps, this shows the recessed ground electrode. 
(B) An SEM cross section of the layered SOI 

structure at the interface between two 
compensation electrodes. The deep V-shaped 
undercut into the oxide layer is clearly visible 
extending 60 /im into the Si02 layer. 

The trap was mounted onto a ceramic chip car- 
rier, with gold ribbon wire connecting the electrodes 
to the chip carrier bond pads. To ensure good rf 
grounding 820 pF capacitors were wire bonded to 
the static electrodes. The chip was mounted inside 
a UHV vacuum system at a pressure of 8.0x10"^ 
Pa. Ytterbium ions were trapped by applying an 
rf voltage Vq = 455(3) V at a frequency Q/27r = 
32.2 MHz via a helical coil resonator with Q = 160. 
All other electrodes were held at ground. Atomic 
ovens containing natural ytterbium were ohmically 
heated producing a neutral flux of atoms parallel 
to the trap surface. Atoms were resonantly ionised 
using a two-colour photoionisation process and laser 
cooled, with ion fluorescence imaged onto an electron 
multiplied CCD array. The ion secular frequencies 



along the trap principle axes were measured to be 
{uja:',ujy',ujz')/'27r = (1.58, 1.47, 3.30) ± 0.01 MHz re- 
spectively, for an ion trapped in the top right corner 
of the array. Using these measurements, numeri- 
cal simulations of the trap predict a trap depth of 
0.42(2) eV. The ion lifetime with cooling light was 
~ 90 minutes, likely limited by background collisions 
and > 5 minutes without cooling light. To maximise 
the ion- ion coupling rate for quantum simulations. 



which scales as i 



the rf voltage can be decreased 



to 150 V, corresponding to secular frequencies of 
co'aj//27r ~ ujyf/27T ~ 0.5 MHz. A 3D contour plot 
of the trap potential for the central 11 trap sites is 
shown in Fig. 3A, a 2D contour plot is projected 
below showing the potential at the ion height. 

A 2D ion lattice can be seen in Fig. 3B, six ions 
are trapped in adjacent lattice sites, with the view- 
able area limited by our imaging system. If required 
we can deterministically expand from simulations 
on regular lattices by introducing defects into the 
lattice structure; these may include both a missing 
ion or the inclusion of an additional ion at a lattice 
site. Such defects result in an irregular lattice which 
when investigating Bose-Hubbard physics, such as 
superfluid-Mott insulator transitions, will result in 
site dependant phonon tunneling rates [4]. Trap- 
ping many ions at lattice sites also allows for the 
simulation of spin models with spin greater than 1/2 
[22]. Figure 3C shows a three ion Coulomb crystal 
trapped in a single lattice site. 

We have also demonstrated shuttling of ions be- 
tween different sites of the lattice. An ion is initially 
trapped at site S, as marked in Fig. 3B. By lowering 
the rf voltage to minimise the potential barrier and 
applying control voltages to the surrounding elec- 
trodes, the ion is shuttled to position A. Figure 4 
shows two CCD images of the ion before and after 
shuttling, also shown are the typical shuttling volt- 
ages. To shuttle the ion back to E the polarity of 
the voltage on Comp 3 and Comp 6 is reversed. 

It is possible to determine the feasibility of quan- 
tum simulations that could be performed using this 
lattice ion trap chip. By considering the ratio P of 
the change in the Coulomb force to the change in the 
restoring force of one ion due to the displacement of 
a neighbouring ion [22, 26], we flnd (3 = 8.32 x 10"^. 
This indicates that ion-ion coupling is dominant be- 
tween nearest neighbour ions only. Combined with 
a commercially available laser with a wavelength of 
355 nm and a power of 0.4 W focused down to a 
30 /im beam waist, a coupling strength J ^ 2830 
s~^ is achievable. It is also important to consider 
sources of decoherence that can have adverse effects 
on the fldelity of the quantum simulation. A sig- 
niflcant source of decoherence stems from the heat- 
ing rate of the ion's motional state estimated to be 
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Figure 3: (A) A 3D contour plot of the lattice 
potential for the central 11 trap sites with an rf 
voltage of 455 V. The equipotential surfaces 
correspond to: red = 0.04 eV, green = 0.4 eV, light 
blue = 0.5 eV and blue = 0.6 eV. A 2D contour 
plot of the potential at the ion height is projected 
below, the equipotential lines are separated by 0.04 
eV. (B) A CCD image of six ions trapped 
simultaneously on the lattice is shown, the 
hexagonal shaped electrodes are also visible due to 
laser scatter. The area imaged was limited by our 

current imaging system and marked by a black 
circle over an SEM image of the trap. (C) A CCD 
image of a single lattice site with a deterministically 
produced three-ion Coulomb crystal. 



1.6 X 10-1^ V^m-^Hz-i [28], which is a function of 
the ion-electrode distance. K^^^ is defined as the 
ratio of the heating rate to the coupling rate and 
is ~ 1.3 for our ion trap lattice [26]. Additionally, 
decoherence can result from the laser induced spon- 
taneous emission rate which again can be compared 
to the coupling rate. The ratio of the coupling rate 
to the laser induced spontaneous emission is defined 
as Lsim [26] and is ~ 12 for this laser. While the 
simulation error would be large using these partic- 



> 1 
o> 

03 

-2 



A 




— RF 

— Comp 3 

— Comp 6 

— Comp 1,2,4,5 



4 6 
Time [ms] 



3D 

300 J 
200 I 

CQ 

100 

< 
10 



Figure 4: Two CCD images are shown before and 
after shuttling an ion from lattice site S to A. The 
shuttling voltage profile for this shuttling operation 
is shown. To shuttle the ion back to A the polarity 
of the voltages on Comp 3 and 6 is reversed. 



ular parameters, amending the hexagon radius and 
hexagon centre separation to 37 /im and 80 /im, re- 
spectively, would result in an error of 0.25 (here we 
assume a laser power of 0.24 W, l]/27r = 30 MHz and 
Vb = 36 V resulting in {ux' ^ujy')/27r ^1.7 MHz, an 
ion-electrode distance of 40 /im, P = 2.8 x 10~^, 
J ^ 291 s~^, Ksim ^ 1-55 (for a heating rate of 
9.4 X 10~^^ V^m~^Hz~^ [[28]^ assuming a cryogenic 
system) and L^^^ ^ 2.1). 

The combination of a scalable microfabrication ar- 
chitecture for trapping 2D lattices of ions and suffi- 
cient coupling strengths between lattice sites opens 
up a new class of quantum simulations. A triangular 
lattice of trapped ions, as demonstrated here, can be 
used for example to explore frustrated XY models 
and 2D spin liquids [5] . With simple modifications to 
the trap geometry alternative lattice configurations 
such as square or honeycomb lattices, of the type re- 
cently described for the investigation of gauge fields 
and the fractional quantum Hall effect [6, 7], can be 
formed. Work to improve the ion coupling rates and 
local ion control will further the suitability of this 
trap for quantum simulation, but represent only in- 
cremental fabrication and trap design steps. 

By allowing high voltages, enabling large ion- 
electrode separations and long ion lifetimes on a 
microfabricated and scalable lattice geometry, our 
device provides an exciting platform for the imple- 
mentation of analogue 2D quantum simulations us- 
ing trapped ions. Not only is this lattice suitable 
for quantum simulations however, other applications 
may include the generation of 2D cluster states for 
ion trap quantum computing. 
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I. METHODS 

Surface flashover measurements were performed 
on test chips consisting of gold coated silicon islands 
surrounded by the exposed handle layer. The chips 
were glued, using conductive silver glue, onto a ce- 
ramic chip carrier and placed in high vacuum at a 
pressure of 6 x 10~^ Pa. Radio-frequency voltage 



was applied by connecting a quarter- wave helical coil 
resonator to the vacuum feedthrough. This ensured 
> 90 % coupling between the rf amplifier and chip 
at a resonant frequency of 28.0 ± 0.5 MHz with Q 
= 210 =b 15. The rf voltage was steadily increased 
while observing the sample. Upon breakdown the 
voltage was measured using a capacitively coupled 
probe. 

Radio-frequency flashover was observed to occur 
~ 20 % lower than static flashover. The exact cause 
for this is not known, however it is likely a result of 
local heating on the chip and outgassing from the 
surface when the rf voltage is applied. To apply suf- 
ficient voltage to trigger fiashover, up to 30 W of rf 
power was applied to the chip. This power would 
be predominately dissipated on chip, which has the 
largest impedance. This would result in significant 
heating of the chip and the release of gases previ- 
ously adsorbed to the surface. 



[1] Schneider, C, Porras, D., and Schaetz, T. Rep. 

Prog. Phys. 75(2), 024401 (2012). 
[2] Blatt, R. and Roos, C. F. Nat. Phys. 8, 277-284 

(2012) . 

[3] Raussendorf, R. and Briegel, H. J. Phys. Rev. Lett. 

86, 5188-5191 (2001). 
[4] Johanning, M., Varon, A. F., and Wunderlich, C. J. 

Phys. B: At. Mol. Opt. Phys. 42, 154009 (2009). 
[5] Schmied, R., Roscilde, T., Murg, V., Porras, D., and 

Cirac, J. I. New J. Phys. 10, 045017 (2008). 
[6] Bermudez, A., Schaetz, T., and Porras, D. New J. 

Phys. 14, 053049 (2012). 
[7] Shi, T. and Cirac, J. 1. Phys. Rev. A 87, 013606 

(2013) . 

[8] Bloch, I., Dalibard, J., and Nascimbene, S. Nat. 

Phys. 8, 267-276 (2012). 
[9] Aspuru-Guzik, A. and Walther, P. Nat. Phys. 8, 
285-291 (2012). 

[10] Peng, X., Zhang, J., Du, J., and Suter, D. Phys. 
Rev. Lett. 103, 140501 (2009). 

[11] Wineland, D.J. and Leibfried, D. Laser Phys. Lett. 
8(3), 175-188 (2011). 

[12] Islam, R., Edwards, E. E., Kim, K., Korenblit, S., 
Noh, C, Carmichael, H., Lin, G. D., Duan, L. M., 
Wang, C. C. J., Freericks, J. K., and Monroe, C. 
Nat. Commun. 2, 377 (2011). 

[13] Friedenauer, A., Schmitz, H., Glueckert, J. T., Por- 
ras, D., and Schaetz, T. Nat. Phys. 4, 757-761 
(2008). 

[14] Gerritsma, R., Lanyon, B. P., Kirchmair, G., 
Zahringer, F., Hempel, C., Casanova, J., Garcia- 
RipoU, J. J., Solano, E., Blatt, R., and Roos, C. F. 
Phys. Rev. Lett. 106, 060503 (2011). 

[15] Gerritsma, R., Kirchmair, G., Zahringer, F., Solano, 
E., Blatt, R, and Roos, C. F. Nature 463, 68-71 
(2010). 

[16] Lanyon, B. P., Hempel, C., Nigg, D., Miiller, M., 



Gerritsma, R., Zahringer, F., Schindler, P., Bar- 
reiro, J. T., Rambach, M., Kirchmair, G., Hennrich, 
M., Zoller, P., Blatt, R., and Roos, C. F. Science 
334(6052), 57-61 (2011). 

[17] Kim, K., Korenblit, S., Islam, R., Edwards, E. E., 
Chang, M.-S., Noh, C, Carmichael, H., Lin, G.-D., 
Duan, L.-M., Wang, C. C. J., Freericks, J. K., and 
Monroe, C. New J. Phys. 13, 105003 (2011). 

[18] Islam, R., Senko, C, Campbell, W. C, Korenblit, 
S., Smith, J., Lee, A., Edwards, E. E., Wang, C. 
C. J., Freericks, J. K., and Monroe, C. arXiv: 1210. 
0142 (2012). 

[19] Korenblit, S., Kafri, D., Campbell, W. C, Islam, 
R., Edwards, E. E., Gong, Z. X., Lin, G. D., Duan, 
L. M., Kim, J., and Monroe, C. New J. Phys. 14, 
095024 (2012). 

[20] Britton, J. W., Sawyer, B. C, Keith, A. C, Wang, 
C. C. J., Freericks, J. K., Uys, H., Biercuk, M. J., 
and Bollinger, J. J. Nature 484, 489-492 (2012). 

[21] Cirac, J. I. and Zoller, P. Nature 404, 579-581 
(2000). 

[22] Porras, D. and Cirac, J. I. Phys. Rev. Lett. 92, 

207901 (2004). 
[23] Chiaverini, J. and Lybarger, W. E. Phys. Rev. A 

77, 022324 (2008). 
[24] Kumph, M., Brownnutt, M., and Blatt, R. New J. 

Phys. 13(073043) (2011). 
[25] Clark, R. J., Lin, T., Brown, K. R., and Chuang, 

1. L. J. Appl. Phys. 105(1), 013114 (2009). 
[26] Siverns, J. D., Weidt, S., Lake, K., Lekitsch, B., 

Hughes, M. D., and Hensinger, W. K. New J. Phys. 

14, 085009 (2012). 
[27] Schmied, R., Wesenberg, J. H., and Leibfried, D. 

Phys. Rev. Lett. 102, 233002 (2009). 
[28] Hughes, M. D., Lekitsch, B., Broersma, J. A., and 

Hensinger, W. K. Contemp. Phys. 52(6), 505-529 

(2011). 



6 



[29] Britton, J., Leibfried, D., Beall, J. A., Blakestad, sen, M. J., Schwab, K., and Monroe, C. Nat. Phys. 

R. B., Wesenberg, J. H., and Wineland, D. J. Appl. 2, 36-39 (2006). 

Phys. Lett. 95(17), 173102 (2009). 
[30] Stick, D., Hensinger, W. K., Olmschenk, S., Mad- 



